Communication through plasma sheaths 



A. 0. Korotkevich/'B A. C. Newell, 2 'Q and V. E. Zakharov 3 - 4 ' ^ 5 '0 
1 Landau Institute for Theoretical Physics RAS, 
2, Kosygin Str., Moscow, 119334, Russian Federation 
2 Department of Mathematics, The University of Arizona, 
617 N. Santa Rita Ave., Tucson, AZ 85721, USA 
3 Department of Mathematics, University of Arizona, 
617 N. Santa Rita Ave., Tucson, AZ 85721, USA 
^Lebedev Physical Institute RAS, 
53, Leninsky Prosp., GSP-1 Moscow, 119991, Russian Federation 
5 Waves and Solitons LLC, 918 W. Windsong Dr., Phoenix, AZ 85045, USA 
(Dated: February 1, 2008) 

We wish to transmit messages to and from a hypersonic vehicle around which a plasma sheath 
has formed. For long distance transmission, the signal carrying these messages must be necessarily 
low frequency, typically 2 GHz, to which the plasma sheath is opaque. The idea is to use the plasma 
properties to make the plasma sheath appear transparent. 



I. INTRODUCTION. 



A. General discussion. 



A vehicle moving through the stratosphere (altitudes 
40km-50km) at hypersonic velocities (8-15 Mach) is cov- 
ered by a plasma sheath. Typically, the plasma density 
n can be as high as 10 18 m -3 with corresponding plasma 
frequency 
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of about 9GHz. In |T]), e is the electron charge — 1.6 x 
1CT 19 C, e = 8.85 x IQ-^CV^m- 1 and M is the elec- 
tron mass 9 x 10 -31 kg. Therefore the plasma is opaque 
to frequencies lower than 9GHz. Direct communication 
through such a plasma to and from the vehicle is im- 
possible because frequencies / suitable for long distance 
propagation through the atmosphere are usually much 
less. For example, the standard frequency used for nav- 
igational satellite systems, including the global position- 
ing system (GPS), are less than 2GHz. For the GPS, 
/ = 1.57542Gffz. 

The challenge is to devise means to maintain contin- 
uous contact with the hypersonic vehicle. When such 
vehicles were principally spacecrafts, a blackout period 
of up to two minutes was acceptable albeit undesirable. 
But when the vehicles are of military origin, it is clear 
that continuous contact is essential for both targeting 
and rapid abort reasons. 

It is a challenge which has drawn many responses. 
They fall into several categories. The first ignores the 
presence of the plasma by using signals with frequencies 
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well above the plasma frequency. The difficulty with this 
method is that such signals are heavily attenuated in and 
scattered by the atmosphere. A second means, which also 
ignores the plasma, is to use low frequency signals in the 
100 M Hz range where wavelengths are large compared to 
the plasma sheath thickness (typically of the order of a 
meter). But such solutions have high cost and low bit 
rates and are not well supported by existing infrastruc- 
ture. A third category of solutions violates the plasma. 
One approach is to remove, by vehicle reshaping, for ex- 
ample, the plasma from certain points on the vehicle at 
which one might place an antenna. Another is to destroy 
it by electrophilic injection or by injecting water drops. 
A third approach is to use powerful magnets to reshape 
the plasma. Such solutions involve a heavy cost in that 
design features necessary for their implementation must 
be built into the vehicle a priori. Nevertheless some are 
feasible and worthy of consideration. For example, it is 
possible to build an antenna into a sharp leading edge 
which would protrude beyond the plasma and survive 
for sufficiently long (it would be eventually destroyed by 
ablation) to cover the flight time. 

The fourth category of solutions, and the one to which 
we are attracted, uses the properties of the plasma itself 
to affect transmission in the same way a judo expert uses 
the strength and motion of an opponent to defeat him. 
One idea is to create new modes of oscillation and prop- 
agation by the introduction of magnetic fields. Indeed, 
for strong enough fields, the Larmour frequency fharmour 
is sufficiently large that the window {fharmour, niax(/^)) 
for which the plasma is opaque is small and transmission 
can be achieved for frequencies below fharmour- But the 
introduction of magnetic fields involves large additional 
weight and new design features. The second idea is much 
more simple. Its aim is to take advantages of nonlinear 
properties of plasma to render it effectively transparent 
to the signal. Communications both to and from the 
vehicle are feasible using basically the same ideas. We 
shall first describe the "to the vehicle" case. Consider 
Figure [1] in which we show schematically the response of 
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FIG. 1: uiL(zr) = lu cos 4>,uil(0) = u>. If the thickness of the 
plasma sheath is equal to L + R — lm, the signal frequency 
/ = 2GHz and the plasma frequency /l — 9GHz then L ~ 
5cm and R ~ 95cm. 



the plasma to an incoming signal with low frequency lo 
from a direction which makes an angle <j> with the nor- 
mal to the vehicle. There are two principal features to 
the response. First, there is a reflection from the layer 
at a point z = z r where the plasma frequency at the 
point LOi,(z r ) is uj cos (f). However, the influence of the 
signal is felt beyond that point, namely at the resonant 
layer z = where wl(0) — Langmuir oscillations are 
excited there which produce large transversal and longi- 
tudinal components of the electric field. The resonant 
layer acts as an antenna. The task is to find a way to 
connect the antenna at the resonant layer at z — to a 
receiver on board the vehicle at z = R. There are several 
possibilities which we have outlined before [H, 0, H| • 

The most practical one, however, is also the most sim- 
ple and first suggested without a detailed numerical sim- 
ulation in We use an onboard source, which we call 
the pump, to generate electromagnetic signals of suffi- 
ciently high frequency lo p (lj p > maxwi(z) + lu) that 

they can propagate through the plasma. There are sev- 
eral candidates for such a source. For example, available 
on the open market is a klystron amplifier which can 
generate ikW of power at frequencies of 12 — \AGHz. 
These high frequency waves have only to travel distances 
of a meter or less. They interact nonlinearly with and 
scatter off the signal wave. Not surprisingly, the largest 
contribution to the scattered wave comes from the non- 
linear interaction of pump wave with the plasma density 
distortion induced by the incoming signal wave at the 
resonant layer. We call the scattered wave a Stokes wave 
because the scattering process is a three wave interaction 
analogous to Raman scattering. The Stokes wave with 
frequency los — lo p —lo carries the information encoded on 
the signal wave back to the vehicle. We will show that, 
whereas much of the scattered Stokes wave propagates 
away from the vehicle, a significant fraction is returned 
to the vehicle. 

What is remarkable is this. The ratio of the power flux 
of the Stokes wave received at the vehicle to the power 
flux contained in the signal wave at the plasma edge can 
be between 0.7 and 2 percent. This means that recep- 



tion of GPS signals may be possible because one simply 
needs an onboard receiver approximately 100 times more 
sensitive than commercially available hand held receivers 
or use sufficiently larger antenna. We shall discuss in the 
conclusion the sensitivity required for a variety of sources. 

Communications from the vehicle requires two power 
sources on the vehicle. One, which we term the Stokes 
wave generator, will also carry the signal. The other 
is the pump wave. Both have carrier frequencies above 
that of the maximum of the plasma frequency. Their 
nonlinear interaction in the plasma produces an oscilla- 
tions of frequency lo = uj p — los ■ Consider Figure O For 
z r < z < R where z r is determined by Wi(z r ) = lo cos 4> 
and (j) is calculated from the differences in propagation 
directions of the pump and Stokes waves, the oscillation 
does not propagate and its strength decays away from the 
vehicle. Nevertheless this oscillation is sufficiently strong 
to act as a power source for a propagating wave in the 
region z < z r where lo cos <f> > lol{z). In the conclusion 
we analyze what power is required in order for the sig- 
nal to be detected by distant receivers. It appears that 
even if we use usual available on the market generators 
communication can be put into practice. 



"Stokes" wave 




Pump wave 



FIG. 2: The concept for communication from the vehicle. 
Although drawn in such a way that the angles of pumping, 
Stokes and signal waves are different, the optimal configura- 
tion is when all angles are the same, i.e. Stokes and pump 
waves are generated in the same direction as the target of the 
desired low frequency signal. 



B. Plan of the paper. 

The plan of the paper is as follows. We begin in Sec- 
tion 2 with a detailed analysis of the two dimensional 
propagation and interaction of a signal wave of frequency 
lo, a pump wave of frequency lu p and a Stokes wave of fre- 
quency us through a plasma with a given density profile 
no(z) where z is the direction normal to the vehicle. The 
key equation is a modification of the well known Ginzburg 
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for the magnetic field amplitude (H(y, z), 0, 0)e" lf2t of an 
oscillation of frequency f2. In the effective electric 
susceptibility is 



s(z,Q) = e Q [l- 
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(ul{z) is the local plasma frequency and v the collision 
frequency). The susceptibility is due to the linear re- 
sponse of the plasma to the electric fields of whichever 
waves are involved. The nonlinear current jnl will be de- 
termined both by the product of the plasma density dis- 
tortion with the linear current and the nonlinear response 
of the electric velocity field due principally to dynamic 
pressure forces. We observe that, for Q ^> maxcji(z), 

z 

the electric susceptibility is approximately £o and the left 
hand side of the nonlinear Ginzburg equation ^ is the 
usual wave operator. 

How do we use ([1])? For the case of communication to 
the vehicle, we use it in two ways. First with j'jvl = 0, 
we determine for Q = w and H(y,z) = H{z)e i{ - u /^ vsiu< t" , 
the field H(z) from which the distortion to the plasma 
produced by the incoming wave is calculated. In this 
instance, H(z) satisfies 



d 2 H 1 de(z,uj)dH 



dz 2 e(z,uj) 
so 



dz 

,2 i 



dz 
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A glance at the third term shows that propagation is im- 
possible for e/eo < sin 2 ^ or, from (O, for wcos0 < 
u>l(z). The importance of the resonance layer where 
e(z,ui) ~ is seen from the denominator in the second 
term. Having solved for H{z) from (H|) we can then cal- 
culate the plasma distortion field 5n(z). Its interaction 
with the pumping wave then produces a nonlinear current 
jn l which gives rise to the Stokes wave. The Stokes wave 
Hs(y, z) and its propagation is calculated by solving @ 
with this jml and appropriate boundary conditions at 
the plasma edge and at the vehicle. Our goal is to de- 
termine Hs(y,z = R). We give the results of both the 
numerical simulation and an analytic estimation. The 
latter takes advantage of the fact, that, for the Stokes 
wave, los S> maxwi(z) and that the principal plasma 

z 

distortion occurs at the resonance layer. 

For communicating from the vehicle, we solve (|4]) with 
the right hand side given by — V x ^tJnl with j NL cal- 
culated from the nonlinear interaction of the pump and 
Stokes waves. Here the goal is to calculate the flux of 



power of the signal wave with frequency u> = u p — los as 
it leaves the plasma edge in the direction of some distant 
receiver. 

In Section 3, we describe the numerical procedure and 
give detailed results of our calculations. 

Finally, in Conclusion, we use our results to calculate 
the powers of both the incoming and outgoing signals at 
their respective receivers. We discuss in addition several 
important considerations: 

• The advantages, particularly in terms of available 
power, of using pulsed signals. 

• The possibility of using GPS sources for incoming 
signals. 

• The challenges involved in making ideas practica- 
ble. 



II. ANALYTICS. 
A. Basic theory. 

We shall study a very idealized situation when the 
plasma sheath is a flat slab. The plasma density is a 
linear function of the horizontal coordinate z 



n (z) = n 



z + L 
R + L' 
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In this geometry the vehicle is the vertical wall placed at 
z = R. The plasma density near the vehicle is no- The 
plasma contacts the vacuum at z — —L, where n — 0. We 
shall study two situations: communication to the vehi- 
cle and communication from the vehicle. In both cases, 
three almost monochromatic electromagnetic waves ex- 
ist in plasma. Two of them have high frequencies lu p 
(pumping wave), cus (Stokes wave). The third one has 
low frequency to, satisfying the condition 



UJ = LOp — LOS- 
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In the "to the vehicle" case w is the circular frequency 
of the incoming signal. In the "from the vehicle" case, u> 
is the circular frequency of the outgoing signal. In both 
these cases, the low-frequency signal plays a key role. 
Because the local plasma frequency at z = is uj, 
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Let us denote also the Langmuir frequency at the vehicle 

as 

2 e 2 «o 
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In a realistic situation /& ~ 9GHz (it corresponds to 
n Q = 10 18 m- 3 ), / ~ 2GHz, R+L = lm, and L ~ 0.05m. 
The wavelength of the incoming signal in the vacuum 
is A = c/f = 0.15m, so that A > L. We point out 
that in the case of low-frequency wave reflection from 
the ionosphere, the situation is the opposite A << L. 

We shall assume that the ions' positions are fixed and 
the plasma is cold (T e ~ 0). The magnetic field has only 
one component H x . The electric field has two compo- 
nents E y , E z . Neither the electric nor magnetic fields 
depend on the ^-coordinate. Maxwell's equations read 

E(0, E y (y, z), E z (y, z)); H(H(y, z),0, 0) 



rv ft 9H 
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V • eqE — e(n — no(z)), j — env. 
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In our geometry, (|14p is one scalar equation. We should 
stress that this is an exact equation. The only challenge 
is the calculation of jn l ■ 

Finally for the magnetic field, one obtain the Ginzburg 
equation 
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For the high frequency pump and Stokes waves e ~ eq. 
Some exact solutions of simplified versions of the homo- 
geneous Ginzburg equation for several important cases 
can be found in Appendix [XJ 

What we are going to do is the following: in subsection 
IIIBI we shall calculate linear responses of the plasma to 
an electromagnetic wave, such as electron velocity, linear 
current and the electron density profile perturbation; the 
calculation of the first nonlinear correction to the linear 
current is done in subsection lll CI analytic estimations for 
"to the vehicle" and "from the vehicle" cases are given 
in subsections III Dl and III El respectively. 



B. Linear responses. 

In order to calculate the nonlinear current we need to 
consider the linear responses of the plasma to the pres- 
ence of an electromagnetic wave. For a field with fre- 
quency f2 
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The power flux in vacuum is 
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In equation (TTB"]) v is the effective friction of the elec- 
tron fluid with the neutral gas, sometimes called the ion 
collision frequency. We take v — 10 8 Hz. 

The current j = jh + jNL- jL is the linear response 
of the plasma on the electric field, jNL is the current 
due to nonlinear effects. For a monochromatic wave of 
frequency f2, Maxwell's equations can be rewritten in the 
following form 
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Using Maxwell equations one can express all responses in 
terms of magnetic field 

fi =ssk(°' (17) 

e 1(8 d \ 

' L = _ M0 2 e(fi) 1 + \ >di H '-dji H )> (18) 

(19) 
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The expression for a distortion Sn of the electron density 
in the plasma n(z) = n$(z) + Sn(y,z,t) can be derived 
from (IS]) and |l8|), 
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C. Nonlinear current. 

The nonlinear current is due to the first nonlinear cor- 
rection to the linear response velocities of electrons and 
a scattering of an electromagnetic wave on the distortion 
of the charge density profile produced by another wave 



Jnl = en (z)vNL + eSnvL- 



(21) 



We introduce the nonlinear velocity wjvl which can be 
found from the following equation 

= vjj x [V x v L ] + i^v L xH--Vv 2 L = -fJv\. 

Here we used a corollary of the Maxwell equations and 
(fT6|) whence to within Oiyjui), 

[Vx^]=-^. 

This means that only the dynamic pressure induced by 
the fields affects the plasma. 

Finally, we have everything for the calculation of the 
first term in right hand side of Ginzburg equation (|15D 
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The detailed expression of equation's (TT5]) right hand side 
can be found in Appendix [Bj 



D. Analytic estimation. "To the vehicle." 



We would like to estimate the ratio 
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S s (z = R) 



of the fluxes of the squared scattered field to the squared 
incoming signal field and express it as a function of pump 
power flux S p measured in Watts per square meter. 

We can make an analytic estimation of the three- wave 
process efficiency. The main contribution comes from the 
vicinity of z = 0. The reason comes from the fact that 
the real part of dielectric susceptibility © for the low 
frequency signal wave has a zero at this point. It means 
that our nonlinear current on the right hand side of the 
Ginzburg equation has a very sharp peak near z = 0. A 
typical plot of right hand side is given in Figure [3] This 
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FIG. 3: The typical right hand side (absolute value) of the 
Ginzburg equation in the "to the vehicle" case. Logarithmic 
scale. One can see that the main contribution comes from the 
region of the point z = 0. 



issue is discussed in more detail in Appendix IB 21 

If we consider a high frequency pumping wave we can 
use the plane wave approximation 

H p (y,z,t) = H p e i{ -^ t+k " v - K f z) . 

The low frequency signal wave can be written 

H (y,t) = H(y,z,t)\ z=0 = H(z)\ z ^e^ t+k y\ 

For the Stokes wave, whose frequency is higher than the 
plasma frequency one can use the following approximate 
Ginzburg equation 



dz 



f + n 2 H s = f s , 



(23) 



where fs is calculated from the curl of the nonlinear cur- 
rent given in pi C|) . To solve, we use the method of 
variation of constants. We find 

H s = de 1KSZ +C 2 e~ iKSZ , 
C[e iKSZ + C" 2 e- iKSZ = 0, 
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e tKsv fs(y)dy. 



One can say that C\ is the amplitude of the Stokes wave 
propagating to the vehicle and C2 is the amplitude of the 
anti-Stokes wave propagating from the vehicle. The main 
contribution to C\(R) arises from the vicinity of z = 0, 
where fs{ z ) is almost singular 
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FIG. 4: Dependence of Ci(R) on the signal incidence angle 



E. Analytic estimation. "From the vehicle." 
Equation ([2]) can be rewritten in the following form 
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(28) 



— tt-Unl)z- 
e dy 

It is not too surprising that that the dominant contri- 
bution to the RHS of l|28p is the first term and arises 
from the neighborhood of z = 0. Again, just as in the 
"to the vehicle" case, the resonant layer acts as a trans- 
mitting antenna which will beam the message contained 
on the Stokes wave to a distant receiver at frequency 
uj = uj p — lus- In Fig. [5]we verify that indeed the domi- 
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After some simple but tedious calculations (see Appendix 
IB2]) one finds 

CUR) ~ 27ri-^— cos(26)smU)H p H*{0). (24) 

where 9 is the pumping incident angle. 

Details of these calculations are given in Appendix lB 21 
The angular dependence of H(0), which we call p{4>), can 
be calculated numerically by solving the homogeneous 
Ginzburg equation. In Fig. [4j we plot the product psm.<f> 
against <f). At the optimal value cf> ~ 0.5, p(cf>) sine/) ~ 1/4, 



d(R) ~ — cos26H p H*(-L). 

I Mc z Sqc 

Using the expression S p = \H p \ 2 /(eoc), one gets 
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For the optimal values of incidence angles (0 = 0, <f> ~ 
0.5), the given plasma parameters and L ~ 0.05m, one 
gets the following maximum value of the efficiency coef- 
ficient 



p s ^ 0.9 x 10" 
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(27) 



This is consistent with what we obtain by direct numer- 
ical simulation. 
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FIG. 5: The full right hand side (absolute value) of the 
Ginzburg equation in the "from the vehicle" case (solid 
line) together with the expression used in the approximation 
(dashed line). Logarithmic scale. Before the vicinity of point 
z — almost no forcing is present. Almost all contributions 
in the vicinity of the resonant point comes from the term used 
in the approximation. In the propagation region (z < 0), the 
approximation slightly underestimates the right hand side. 

nant contribution comes from the first term on the RHS 
of P5|) and from the neighborhood of z = 0. Hence we 
can get simple equation for a very good approximation 
to the approximate particular solution of (|28|l . namely 



dH 
~dz~ 



Unl)v 

The general solution is the following 



H = dfa{z) + C 2 Mz) + J {jN L ) y dz, 





(29) 



(30) 



where 4>i{z) and 4>2{z) are solutions of the homoge- 
neous part of equation (|28|) , and <f>i(z) is bounded as 
z — > R 1, 02 (z) is unbounded (exponentially) at the 
vehicle. Thus C 2 — 0. See Appendix lAl for a discussion 
of solutions to the homogeneous Ginzburg equation. 
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Using the boundary condition on the edge of the 
plasma (z = —L) 



dH 

dz 



(-L) = -iKoff(-L), 



where kq = — cos 6 is the z-component of wavevector of 

the outgoing low frequency signal wave, and Jnl{—L) = 
0, one finds 
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Finally for the magnetic field at z = — L we find 

-L 

H{—L) ~ — T t [ ^ L [ , r , [(jNL)ydz. ,:!2 ) 
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The function (jjvl)j/ oscillates with z with wavenumbcr 
Kp — k-s- The lower the wavenumber the more will be the 
contribution in the integral. This gives us a very sim- 
ple optimal strategy for the choice of pump and Stokes 
wave directions. We should radiate both the Stokes and 
pumping waves in the desired direction of the signal wave 
propagation. In this case we also have an exact compat- 
ibility with the boundary conditions at z = — L. 

If we consider the expression for (Jnl)v given in Ap- 
pendix [Bl we can see that in the case ujq -C lus,uj p the 
first term (|B5|) is the major one in the vicinity of resonant 
layer. The resonant layer works like radiating antenna. 

Using the simplified nonlinear current expression and 
considering the pumping and Stokes waves as plane waves 
one finds 



H(-L) ~ - 



eu} L sin 6 1 

„3„._,.. ~ H p H S x 



(f)' x {-L) + iK (t>x{-L) 



2M£ c 3 wsw p A 

qiA cos <{> 

j4cos< 



(33) 



Where A = Laj a /c. 

Using the solutions of the approximate homogeneous 
equations (|A8[> . we can estimate fli(z)/(l>i(z)\z=-L — 
1/L. Thus for kqL = Acosd <C 1, one finds 



H(-L) 
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For the power density we have 



32 \Mc 2 J 



e c \uj s u! p 



sin 2 (j)SsSp. (34) 



This result is quite clear from physical point of view. The 
larger is, the longer is the distance over which the signal 
wave is generated in the plasma. 



In our simulations, A ~ 2.1 and in this case we cannot 
use the simplified expression given above. Instead we 
find, 
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1 + Cder COS 2 . 

Here we introduced the coefficient Cder = ( K o0i/Vi) 2 the 
value of which we obtain from our numerics. 
Finally, we find 

sin(^4 cos 6) 
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i. r,.,,, 2 , o sin(Acos0) 



^4 cos< 



f 1 — cos(A cos < 
A 2 cos 2 



1 



1 + Cder COS 2 



S S S P . (37) 



The subscripts refer to the frequencies of the onboard 
pump waves. Again, we find the magnitude and angular 
dependence to be consistent with our numerical results. 



III. NUMERICAL PROCEDURES AND 
SIMULATIONS. 

The equation we solve numerically in all cases is the 
Ginzburg equation (|15p including all terms on its right 
hand side. The boundary conditions are given at z = 
L\ = —L — {L + R), in the vacuum beyond the plasma 
edge and at z = R, the vehicle. 

To solve this equation we use a "sweep" -method de- 
scribed in detail in Appendix [Cj The method was in- 
vented simultaneously in several places for work on classi- 
fied topics in the middle of the last century. In the Soviet 
Union, it was introduced by a group of L. D. Landau (in- 
formation from I.M. Khalatnikov) (the first publication 
[5| appeared several years later due to obvious reasons) 
and was developed to its modern form in Q. 

As the first step in the "to the vehicle" case we have 
to find the profile of the incoming magnetic field in the 
plasma. We used an incident angle 6 = 0.5. It will 
be shown later that this angle is an optimal value but 
it is good for an initial evaluation of the possibility of 
communication. We consider the incoming signal as a 
monochromatic plane wave of a given frequency /o = 
2GHz and amplitude Ho. The current is equal to zero. 
In this case, the boundary conditions are 



dH 

-L\, — 1- ikqH = 2iK H Ql 

oz 



R. 



dH 

dz 



0. 



(38) 
(39) 
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The resulting profile of the magnetic field is shown in 
Fig. El The profile of E z (z) is shown in Fig. At 
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FIG. 6: Incoming signal magnetic field profile. 
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FIG. 8: Magnetic field profile of the Stokes wave. Pumping 
frequency 12GHz. 

1.5x10 7 



1.0x10 



5.0x10"° - 



0.0x10 



-5.0x10 



-1.0x10" - 



-1.5x10 




-1.2 -1 -0.8 -0.6 -0.4 -0.2 0.2 0.4 0.6 0.8 1 



FIG. 7: Incoming signal electric field (z-component) profile. 



FIG. 9: Magnetic field profile of the Stokes wave. Pumping 
frequency 18GHz. 



the next stage, we consider an incident low frequency 
magnetic field profile as a source of distortion of the 
plasma density profile and take into account currents due 
to the presence of a pump wave. The pumping wave angle 
6 = 0.0. Our goal is to calculate the scattered field Hs 
with frequency u>s = u p — u>. In this case, the boundary 
conditions are 



= -u 



dHs 

dz 



ik s H = 0, 



OZ 



(40) 
(41) 



The profiles of the magnetic fields Hs for two different 
pumping frequencies are shown in Figures [8] and [9l We 
note that the resonant layer z = acts as if it were a 
source. 

In the "from the vehicle" case we calculate the mag- 
netic field of the low frequency wave generated by plane 
pump and Stokes waves. Following the optimal strategy 
in this case, described in the analytic part of the paper, 



we take all angles equal to each other <j> — 
this case, the boundary conditions are 

dH 

z = -Li, — h ikqH = 0, 

oz 

z = R, H = 0. 



7r/4. In 

(42) 
(43) 



Here H(z) is the magnetic field of the signal wave with 
frequency w = w p - wg. The boundary condition at z = 
i?, H = 0, gives us the worst of all cases by definition. 

The low frequency magnetic fields for two different 
pumping frequencies are shown in Figs. [TUland llll 

We tested the robustness of the code by allowing for 
both finite and zero conductivity of the vehicle surface in 
the "to the vehicle case". During the simulation in the 
"from the vehicle" case we also redid the simulation with 
the derivative of the magnetic field at the vehicle equal 
to zero. In all the cases, the influences of the differing 
boundary conditions were negligible. 

In the "to the vehicle" case, it is convenient to intro- 
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of the output signal flux to the product of the pump and 
Stokes fluxes and express it as a function of the optimal 
angle. 
We found 



ojp = 2tt * 12GHz, max(/i) ~ 1.8 x 1(T 16 
up = 2tt * 18GHz, max(/i) ~ 3.0 x lCT 17 



1 

lWm~ 2 ' 
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In order to investigate the dependence of the result on 
the angles tf>,9 p ,6g, we calculated fi for various different 
choices. The results are shown in Figs. n^lfTTl 



FIG. 10: Generated low frequency magnetic field. Pumping 
frequency 12GHz. 
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FIG. 12: Dependence of power conversion efficiency coefficient 
fis on angles. "To the vehicle". Pumping frequency 12GHz. 



FIG. 11: Generated low frequency magnetic field. Pumping 
frequency 18GHz. 



duce the function fig a s the ratio 

S s (z = R) 
Ms = ~ 

JO 

of the scattered field flux to the incoming signal flux and 
express it as a function of pump flux S p measured in 
Watts per square meter. We found 



v-s 



lu p = 2tt* 12GHz, max(/j s ) ~ 2.2 x 10" 



2?r * 18&ffz, max(^s) ~ 0.63 x 10 



lWm 2 ' 



-li 



S n 



lWm- 2 ' 



These results are in a good agreement with the analytic 
estimation (|27[) . Any difference is due to the fact that the 
pumping frequency is not sufficiently high to neglect the 
plasma frequency. The reason we used these frequencies 
and not much higher ones was that they are available on 
standard microwave equipment and devices. 

In the "from the vehicle" case, we calculate the ratio 

_ S ou t(z — —L) 
SgS p 




FIG. 13: Dependence of power conversion efficiency coefficient 
fis on angles. "To the vehicle". Pumping frequency 18GHz. 

As one can see, in the "to the vehicle" case we have 
a very good agreement between the analytically esti- 
mated angular dependence (|26[) and the numerical re- 
sults. Namely, we have a maximum at pumping angles 
close to 9 = and the efficiency coefficient fig goes to 
zero at the vicinity of 9 — 7r/4 in a agreement with the 
cos(2#) dependence. So we can formulate a simple rule: 
in order to get the best possible performance, send the 
pump wave in a direction perpendicular to the plasma 
edge surface. 
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FIG. 14: Dependence of the power conversion efficiency coef- 
ficient fis on several pumping angles, in the "to the vehicle" 
case. Pumping frequency 12GHz. 
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FIG. 16: Dependence of the power conversion efficiency co- 
efficient fi on optimal angle. "From the vehicle". Pumping 
frequency 12GHz. 
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FIG. 15: Dependence the of power conversion efficiency coef- FIG. 17: Dependence of power conversion efficiency coefficient 
ficient fis on several pumping angles, in the "to the vehicle" fi on optimal angle. "From the vehicle". Pumping frequency 



case. Pumping frequency 18GHz. 



18GHz. 



In the "from the vehicle" case, the situation is even 
simpler. As it was shown in Section III El the power con- 
version is optimal if we radiate both the pump and Stokes 
waves in the direction of the desired signal wave propa- 
gation. The estimated angular dependence ([35]) can be 
fitted with good accuracy to the numerical results using 
only one tuning coefficient Cder- It is shown that this 
coefficient weakly depends on the pumping frequency. 

IV. CONCLUSION AND DISCUSSION. 

Let us now discuss the practical usage of this approach 
for receiving at and transmitting from the vehicle. For 
the "to the vehicle" case we consider the problem of re- 
ceiving even GPS signals. Let us estimate the resulting 
attenuation coefficient. Given a pump waveguide aper- 
ture of 3cm x 3cto and a pump power of 3kW, this 
gives S p = 3.3 x 10 6 Wm~ 2 . One can use the pulse 
regime. In this case, even for pulses 10~ 3 s long, every 
pulse still contains more than 10° periods of the low fre- 



quency signal and we can get much higher power flux 
gpuise = 3 3 x 10 9iy m -2. i t g i V e S us the attenuation 

coefficients fisS^ ulse 

li s S p p ulse ~ 0.73 x 10~ 2 , uj p = 2tt * \2GHz. 
^ sS pnise _ 2.1 x 10~ 2 , uj p = 2tt * 18GHz. 

The usual level of a GPS signal at the Earth surface 
is about — 127.5dBm (1 Decibel per milliwatt is equal 
to IdBm = 101og 10 (P/lmI^)). Indoors, one must use 
high sensitivity GPS receivers. Many general purpose 
chipsets have been available for several years. Presently, 
the market offers sensitivities —157.5dBm (for example 
7]). Using the definition of dBm one can see, that it 
is possible to receive a signal with an attenuation about 
10~ 3 . Also it is possible to use a much bigger antenna 
on the vehicle than in the case of a handheld device. 
In this case, it is even possible to receive a signal using 
the continuous rather than pulsed regime for a klystron 
pump. So even at the angles far from optimal, one can 
receive GPS signals. Further, we used maximum value 
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FIG. 18: Schematic plot of beam diffraction. 



of the plasma thickness. If the plasma sheath is thinner, 
the angular dependence is broader. 

Available on the open market are klystron 
amplifiers with the following characteristics If: 



Model 


Frequency (GHz) 


Power (kWatt) 


Mass (kg) 


LD4595 


14.0-14.5 


3 


40 


LD7126 


17.3-18.4 


2 


27 



In the "from the vehicle" case, because of sensitive land 
based receivers, all we need is to have a reasonable signal. 
Let us estimate an incoming power on some land based 
antenna. First of all, for any real antenna we have to take 
into account the decrease of a signal due to diffraction 
broadening. If the diameter of the land-based antenna 
(Figure [18]) is Dq, the diameter of the signal flux after 
some long distance I will be 



D(l) 



2Da 



(44) 



It means that if we have power flux at an antenna Sa, 
the power flux at the edge of the plasma after a distance 
I will be 



(45) 



For example, for an antenna of the diameter equal to 5m, 
after 100km 

So ^ 1.1 x 10" 5 SU. 

Now one can calculate the sensitivity of the receiver 
needed. Let us suppose that the signal beam outgo- 
ing from the vehicle has diameter Dq = lm, signal 
frequency / = 2GHz and corresponding wave length 
A = 1.5 x 10 _1 m, the land based antenna has a diame- 
ter Dlb = 5m and is situated at a distance I = lOOfcm. 
Using the previous results for diffraction, the pumping 
klystrons' powers from the table above and the expres- 
sion S ou t = /iSpSs, one can get for the power on the land 
based receiver 



S 



LB 



(2D% 
V IX 



= 1.8 x lO" 8 ^ 



(46) 



We now list for two different frequencies, the correspond- 
ing powers in Watts at the receiving antenna. 

uj p = 2ir * 12GHz, 

P A ~ 1.8 x 10" 8 * 1.8 x 10~ 16 * 9 x 10 6 VFm~ 2 * 25m 2 

~ 0.73 x KT 15 W; 

ujp = 2tt * 18GHz, 

P A ^ 1.8 x 10~ 8 * 3.0 x 10~ 17 * 4 x 10 12 TUm~ 2 * 25m 2 

~ 0.54 x 10- 17 W. 

The GPS receiver mentioned above has a sensitivity 
about — 16QdBm ~ 10~ 19 VU Even with such a modest 
size of the antenna and ordinary klystrons one can receive 
the signal at almost any angle. 

As a final remark one can conclude that proposed 
method for communication with and from the supersonic 
vehicle is realistic even using standard devices available 
on the open market. 
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APPENDIX A: ANALYTIC SOLUTIONS OF THE 
GINZBURG EQUATION IN SOME SPECIAL 
CASES. 

By neglecting ]nl, we obtain the linear Ginzburg 
equation. It takes an especially simple form if O = u>, 
v/oj = and H ~ e' ky . In this case e/eq = ~z/L and 
equation ([15]) is 



d 2 g 
dl 2 " 



1 dH 

z dz 



Here 



A= -^(L + R) 



z 

A3 

1/3 



H = 0. 



(Al) 



1/3 



, (A2) 
2vr x 9GHz, 



and A is another length. In our case ujl 
R + L = lm and A = 0.03m ~ L. 

One can introduce the dimensionless variable £ = zj A 
Then equation (|A1[) simplifies to, 



d 2 H 1 dH 



- (£ + a 2 )H = 0. 



(A3) 



d£ 2 £ d£ 

Here a 2 — A 2 k 2 is a dimensionless constant. 

Equation (|A1[) has two linearly independent solutions 
<t>i, 4>2- We assume 



0, 



oo, at z — > oo. 
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The Wronskian of these solutions is proportional to e/eq. 
We can put 

W = {fa, fa} = M2 ~ Mi = -j. (A4) 
It means that 

W\ Z= _ L = 1. 

Equation (|A3[) cannot be solved in terms of any known 
special functions. In the "outer" area |£| ^> a 2 it reduces 
to the form 



d 2 H 1 dH 



•£ff = 0. 



(A5) 



One can check that equation (|A5|) can be solved in terms 
of the Airy functions Ai and Bi. Namely, 



1=ai Ai'(O^^ 1/4 e - 2/3 « 3/2 , 
b 2 = b 1 Bi , (0~^=e /4 e 2/3e/2 , ate -00. 

\/7T 



(A6) 



Equation (|A8|) can be solved in terms of Bessel functions 
[4| . Two linearly independant solutions of equation (|A8p 
-01, t/>2 behave in neighborhood of £ = as follows 



^2 



= 1 



-e 4 + 



(A9) 



From (|A4p one gets 



Both solutions, which are some linear combinations of 
ipi, ip2 are bounded. Thus the magnetic field has no 
singularity at z = 0. 



APPENDIX B: RIGHT HAND SIDE OF THE 
GINZBURG EQUATION. 

1. General case. 

Consider the Ginzburg equation for a wave 
H 3 (y,z,t)^H 3 (z)e^ t+ik ^, 



ttA 2 



d 2 H 1 dH 



de c 



a 2 H = 0. 



(A7) 



In the "inner" area |£| -C a 2 , the equation (|A3|) is 
reduced to the form 



(A8) 



and calculate the right hand side of (| 1 5[) in terms if the 
fields Hi(y, z, t), H2(y, z, t). In the "to the vehicle" case, 
Hi will represent the pump wave, Hi the signal wave and 
H 3 the Stokes wave. In the "from the vehicle" case, H 3 
will be the signal and Hi and H2 the pump and signal 
carrying Stokes waves respectively. In all cases uj 3 = 
lu 2 — lux, k 3 = k 2 — k\ . We find 



V x j NL (Hi,H 2 ,ki : k 2 ,k3,uji,uj2,uj 3 ) 

J X 

1 dH* dH 2 

2M 2 u 3 (1 + iv/wa) \sl (1 - w/oji) uj\e 2 (1 + \v/oj 2 ) to 2 dz dz 

k\k 2 



e 6 n' Q (z)k 3 



e* (1 — uv/Vl) u 2 e 2 (1 + iv/oj 2 ) 



fefc 2 



d ( n (z) 



M 2 \ (1 — iv/uji) uj\e 2 (1 + \v/w 2 ) ui 2 dz \ e\ 



Jvn k\ 



d ( n (z) 



(1 - \v/u)i) Ld 2 el (1 + w/uj 2 ) J\ dz \ e 2 
1 dH 2 



H'\H 2 



ki d ( d ( n {z) 

ti 1 



M 2 \uj 2 £ 2 (1 + iv/uj 2 ) dz (1 - iu/ui) dz \ 1 dz \ e\ 

1 dH{ k 2 d f d fn (z) 

H 2 



oj\e*i (1 — \v/u>i) dz 0J2 (1 + \vjw 2 ) dz \ * dz \ e 2 
Moe 3 / k 2 d ( n Q (z)\ ki d ( n (z) 



M 2 \u^(l + iiy/u 2 ) dz V £2 



jf (1 - iv/wx) dz \ e* 



H\H 2 . 



Using formulae (121 



(Jnl(Hx 



H 2 ,k 1 ,k 2 ,k 3 ,uj 1 ,uj 2 



' y 



(Bl) 



(B2) 

(B3) 
(B4) 
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e 3 n (z)k 3 



1 



dHl dH 2 



2M 2 u 3 (1 + iv/uoz) \e* (1 - iv/uJx) u 2 e 2 (1 + if/o^) <9z 9z 



(1 - iu/wi) J\e 2 (1 + WW2) < 



;H*H 2 



1 



<9ff 2 fci 



„ 9 / n (z) 



M 2 \u>2^2 (1 + if/wa) <9z (1 - iz//wi) 1 <9z V e* 



+ 



H 2 — 



d ( n (z) 



)\e\ (1 — \v/uji) dz (\ + iv / uj 2 ) dz \ e 2 



(B5) 



(B6) 



2. Approximate right hand side. "To the vehicle" case. 

In the "to the vehicle" case, the main contribution comes from the terms containing poles 

V x j NL (H,H pi k,k p ,k s ,uJ,LL> p ,ujs) ^ 

J X 

d ( n Q (z)* 



M 2 (1 — iv/ojp) Lu 3 s p (1 + iv/LO p ) to 2 dz 



H*H, 



1 



9 f H * d ( n (z) 



M 2 uj 2 Sp (1 + iv/ujp) dz uj 3 (1 — iv/uj) dz \ dz \ 



k d f no(z) 



M 2 lu 3 (1 - w/w) dz 



H*H, 



v 



r 



(B7) 



(B8) 



Assume that the high frequency pumping wave remains 
undisturbed. Then 



H p (y,z,t)=H p e^ k "y- K " z - 



we find 



and 



fs(z) 



C^R) 



2lK S 



V x j NL 



fs(y)e- iKsy dy 



2lK S 



fs(y)e- iKsy dy. 



After several integrations by parts in the second term of 



V x j NL 
ks)k p 



taking into account k$ = k p — k, one finds 



-ik e 3 H, D H* 



2k s M 2 e^Iuj 3 



(kp(k p - k s ) + (n p + (B9) 



+ OO 



<~ 1 j dz 

-00 



d_ ( Mz)_\ p -i(K S+Kp ) d ^ 



For most pumping angles, and using the fact that uj p ^> 
to, one can substitute u>s — and consider incidence 



angles of pumping and Stokes wave to be close in absolute 
value. Following Fig. [1] the pumping incidence angle is 9 
and low frequency signal incidence angle is (f>. 

Ci (R) ~ * 3H * - 4 cos(2tf) (BIO) 



2k s M 2 e uj 2 uj 3 c 2 



d ( n (z) 



dz V £ 



H* e -^ K3+K ^dz. 



Using integration by parts once more one can get 



Ci{R) 



(k p + Ks) e 3 H } 
' 2 e c 

'n {z) 



2ks M 2 soc 3 lo- 



■ cos(26») sin(0) x (Bll) 



Calculating this integral by residues and taking into ac- 
count 

L uj 2 
n {z) ~ n (z = 0) = n Q - ■ - = n — , 



£0 

e 



we finally get 



d(R) ~ 2iri 



e 3 ngL 



M 2 e 2 c 3 J 1 L 



L + R u) A L 

L A - v t 

z + to UJ 



cos(2<2) sm((f>)H p H*(0) 



27ri-^ — cos(26») sm(6)H v H*(Q). 

Mc z £qC 
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APPENDIX C: NUMERICAL METHOD. 

Here we briefly present formulae for the "sweep" 
method in a very general way following the approach 
given in Q. 



The result is a tridiagonal matrix ((N + 1) x (N + 1)) 
equation for a, b, c. 



2. "Sweep"-method. 



1. Reformulation of a problem on a grid. 



Consider the ODE 

p(x)pL+q(x)%-+r(x)y = f{x), (CI) 
aar dx 

in the region < x < L with boundary conditions 



dy 

a dx + " y \ x =° = 7 ' 

ai dx " iy \ x=L = 71 ' 



(C2) 



We use for (|Clj) a second order finite difference scheme 
on a grid of (N + l)-nodes (yo = y(0), j/jv = y{L)) with 
constant step h 



Vn+i - 2y n + y„-i y n+ i - y n -i _ 

Pn j n "T Qn n , T" ^*n?/n In 



h 2 



2h 



(C3) 



This equation is only valid for inner nodes of the grid. 
The boundary conditions take the form 



Vi-Vo . a 
a 7 H PVo = 1, 



yN ~ VN-l . „ 

ai r 1" PWn = 7i- 



(C4) 



We can rewrite equation (|C3|) as 



a n y n -i - b n y n + c n y n+ i = d n , (C5) 

Pn Qn Pn Qn 

In the same way for (|C4[) . one finds 



&o2/o + co2/i = ^o, (C6) 
aNVN-i - b N y N = d N , (C7) 

(C8) 



u a a a a 
bo = - - p,Co = -,d = 7 

h h 



a N = - — ,o N = -— Pi,a N = 7i 

h h 



The solution of the linear system of equations with 
tridiagonal matrix is well described in numerous sources 
(for instance [HI). It can be shown that one can find a 
solution in the following form 

2/n-l = PnVn + Qn- (C9) 

From the left boundary, we have from (|C6[) . that 

co . d 
yo = t-Vi + T- 
bo b 



In this case 



p _ co n _ rfo 
oo b 



(CIO) 



Next we derive a recurrence relation for P„ and Q n . After 
substituting flCS) in jC5|, we find 



Un 



-Vn+1 



J n ^tf n 



Then, comparing with (|C9[) . we see that 



> Qn+1 



&71 CLn Pn 



(Cll) 



Using the initial values (|C10p and the recurring relations 
(|C11I) , one can get all P n , Q n coefficients up to n = N 
("direct sweep" from left to right). 

Than we use second (right) boundary condition (N-th 
equation) 

aNVN-i - b N y N = d N and y N _i = P N y N + Q N . 
Immediately one finds 

dN — &nQn 



yN 



o-nPn — bjv 



(C12) 



Finally, performing a recurrent "backward sweep" (from 
right to left), using the already known P n ,Q n , "sweep"- 
relations (|(J9[) and the initial condition (|C12[) . we get 
values for all y n . 
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